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I.  INTRODUCTION 


The  University  of  Denver's  part  in  the  Office  of  Naval  Research 
Project  on  Aircraft  Crashworthiness  has  been  directed  at  head  and  neck 
motion  during  crashes.  Specifically,  it  has  been  concerned  with  possible 
constraint  systems  to  protect  the  head  of  an  aircraft  occupant  against 
violent  rotation  about  the  torso  and  subsequent  loss  of  consciousness 
during  emergency  landings  and  survivable  crashes. 

The  restraint  concept  that  has  been  tested  is  an  inflatable 
collar  to  limit  head  motion.  Thp  collar  would  fit  under  the  chin  and 
cushion  the  impact  of  the  chin  on  the  thorax  during  hyperflexion  of  the 
neck.  During  the  first  year  of  the  contract,  an  airbag  was  tested 
using  a  single-degree-of-f reedom  mechanical  system.  An  anthropometric 
dummy  head  was  mounted  on  a  rigid  bar  that  simulated  the  neck.  This 
assembly  was  mounted  on  a  died  that  was  accelerated  by  compressed  air. 

The  results  of  tests  on  this  mechanical  system  showed  that  the 
airbag  is  a  feasible  method  of  limiting  head  motion  during  flexion. 

The  second  part  of  the  contract  has  been  devoted  to  further  explo¬ 
ration  of  tills  concept.  A  mechanical  model  has  been  built  to  better 
simulate  the  human  neck.  This  model  allows  stretching  of  the  neck 
and  relative  angular  motion  between  the  head  and  neck  similar  to  head 
rotation  about  the  occipital  condyles.  This  more  realistic  model 
allows  more  confidence  in  the  experimental  results  evaluating  the 
effectiveness  of  the  airbag. 

In  addition  to  the  experimental  work,  theoretical  work  has  been 
devoted  to  modelling  the  motion  of  the  head  and  neck  during  accelera¬ 
tion.  The  acceleration  levels  of  this  motion  must  be  kept  at  low 
levels  in  order  to  survive  a  crash  without  serious  injury.  Part  of  the 
resistance  to  this  motion  comes  from  the  muscles  of  the  neck  themselves. 

The  muscle  tension  during  the  normal  range  of  motion  affects  the  response. 
After  the  rotation  exceeds  the  voluntary  response  limits,  additional 
torque  is  produced  by  attempting  to  further  stretch  muscles  and  ligaments. 

The  main  theoretical  problem  is  arriving  at  a  tractable  set  of  consti¬ 
tutive  relations  to  compute  neck  moments  as  a  function  of  neck  displace¬ 
ments  and  velocities.  Many  different  representations  have  been  used. 

The  simplest  are  linear  springs  and  dashpots  connected  in  parallel 
or  in  series.  More  sophisticated  models  use  nonlinear  springs,  stops, 
and  viscoelastic  elements.  The  number  of  degrees  of  freedom  varies 
from  one  or  two,  to  models  allowing  three  degrees  of  freedom  for  each 
cervical  vertebrae. 

Muscle  contraction  is  varied  by  the  central  nervous  system.  The 
total  tension  in  the  muscles  depends  on  the  number  of  fibers  excited. 

The  force  does  not  depend  on  the  length  of  the  muscle  during  the  normal 
range  of  motion.  This  independence  of  force  on  length  during  initial 
phases  of  the  motion  indicates  that  linear  springs  are  not  satisfactory 
elements  for  the  mathematical  modelling  of  the  effect  of  muscle  tension. 
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The  recent  literature  has  used  viscoelastic  elements  to  simulate  this 
type  of  muscle  force.  These  elements  contain  viscous  damping  terms. 

A  different  approach  has  been  used  in  the  present  study.  The 
voluntary  resistance  cf  the  neck  muscles  to  motion  is  treated  as  a 
constant  torque  during  the  normal  range  of  rotation.  The  constant 
torque  resisting  the  motion  is  analoguous  mathematically  to  solid 
friction  (or  Coulomb  damping)  rather  than  viscous  damping. 

Nonlinear  hard  springs  are  used  in  the  mathematical  model  to 
simulate  the  increased  torque  when  the  head  and  neck  rotations  exceed 
the  voluntary  limits. 

A  second  feature  of  the  theoretical  part  of  the  study  is  an  attempt 
to  incorporate  a  systematic  approach  for  estimating  the  parameters  in 
the  new  constitutive  relations.  Given  the  motion  of  the  head  and  neck 
from  experimental  data,  the  problem  is  to  deduce  the  equations  of  motion. 
This  estimation  of  parameters  is  treated  by  finding  the  effect  of  vary¬ 
ing  each  parameter  by  Newton's  method  (quasilinearization)  and  applying 
a  least-squares  fit  of  the  experimental  data  to  compute  the  variation 
in  each  parameter.  This  method  of  estimation  theory  has  been  used  for 
process  dynamics  in  chemical  engineering,  and  for  other  applications, 
but  we  have  not  seen  it  used  in  the  field  of  biomechanics  of  human 
motion. 

The  computer  program  to  predict  the  response  for  given  initial 
conditions  and  given  parameters  was  straight  forward  to  write  and  soon 
gave  numerical  results. 

The  estimation  theory  is  an  iterative  procedure  because  of  the  non¬ 
linear  nature  of  the  differential  equations  of  motion.  This  part  of  the 
numerical  computations  showed  poor  convergence.  The  least-squares 
subroutine  to  calculate  changes  in  the  assumed  parameters  in  the  consti¬ 
tutive  equations  contains  a  set  of  linear  algebraic  equations  that  must 
be  solved  to  find  the  corrections  to  the  current  values  of  the  constants 
appearing  in  the  constitutive  equation.  The  determinant  of  this  matrix 
was  nearly  singular.  In  addition,  several  of  the  submatrices  were  also 
nearly  singular. 

The  slow  convergence  of  the  iteration  due  to  linear  dependence 
is  discussed  at  length  in  the  body  of  the  report.  Numerical  solutions 
were  finally  achieved  by  varying  one  parameter  at  a  time. 

The  lack  of  sensitivity  to  parameter  variations  is  of  more  than 
academic  interest.  There  are  at  least  two  physical  explanations  of  the 
linear  dependence  that  are  of  major  significance.  First  of  all,  the 
lack  of  convergence  can  merely  mean  that  there  is  not  a  unique  set  of 
parameters  that  minimizes  the  sum  of  the  squares  of  the  differences 
between  theory  and  experiment.  This  implies  that  if  there  is  not  a 
unique  set  of  parameters  for  a  given  theory,  then  more  than  one  theory 
can  be  used  to  explain  the  same  data.  This  line  of  reasoning  leads 
back  to  the  principle  that  one  experiment  can  disprove  a  theory  and  any 
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number  of  experiments  cannot  "prove"  a  theory. 

A  second  explanation  of  the  linear  dependence  is  related  to  the 
first.  The  least-squares  analysis  was  based  on  computed  and  measured 
angular  rotations  of  the  head  and  neck.  The  rotations  are  calculated 
by  integrating  the  equations  of  motion  for  the  angular  acceleration.  An 
integral  of  a  function  is  not  as  sensitive  to  a  change  in  a  parameter 
as  a  function  itself.  A  change  in  momentum  depends  on  the  impulse  and 
different  forces  can  be  integrated  with  time  to  give  the  same  impulse. 

The  logical  step  would  be  to  base  the  least-squares  analysis  on 
measured  angular  accelerations.  However,  practical  problems  presented 
themselves  in  following  this  approach.  A  complete  set  of  plotted 
experimental  data  for  live  subjects  was  available  for  sled  accelerations, 
angular  accelerations  of  the  head  and  neck,  and  head  accelerations  plus 
the  velocities  and  the  displacements  associated  with  these  accelerations. 

However,  the  linear  acceleration  of  the  base  of  the  neck  at 
vertebrae  T-l  was  not  plotted  and  this  acceleration  is  required  for 
input  data  to  the  theorv.  This  T-l  acceleration  appears  to  differ 
markedly  from  the  sled  acceleration  because  of  slack  in  the  shoulder 
harness,  elasticity  of  the  harness  and  of  the  thorax  and  shoulders  of 
the  human  subject. 

The  results  in  this  report  are  computed  using  the  sled  accelera¬ 
tion  in  lieu  of  the  T-l  acceleration.  This  was  sufficient  for  computing 
theoretical  angular  deflections,  but  the  correct  set  of  acceleration  data 
is  necessary  before  proceeding  further  with  the  least-squares  analysis 
of  the  acceleration. 
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II.  CONCLUSIONS  AND  RECOMMENDATIONS 


Before  going  into  the  details  of  the  work  on  the  second  phase  of 
the  contract,  it  might  be  useful  to  outline  the  results  obtained  thus 
far,  how  they  fit  into  the  current  literature  on  head  and  neck  motion, 
and  the  directions  that  future  work  should  take  in  this  field. 

1.  The  mechanical  neck  developed  under  the  contract  has  some 
novel  features  that  add  realism  to  its  function  as  a  means  of 
evaluating  proposed  head  restraint.  The  use  of  shock  cords 
to  simulate  muscle  tension  allows  adjusting  this  variable 
without  major  modification  of  the  model.  A  catch  has  been 
designed  to  oppose  rebound  from  the  fully  flexed  position. 
This  prevents  the  oscillation  observed  in  the  simple  one- 
degree-of-f reedom  model  used  on  the  previous  evaluation 

of  an  airbag  and  is  closer  to  experimental  observations 
on  human  subjects. 

2.  Tne  airbag  acts  to  relieve  the  neck  forces  due  to  flexion. 
This  is  achieved  at  the  price  of  inducing  forces  on  the  jaw 
and  sternum.  Data  are  scarce  for  a  trade-off  study  between 
neck  injuries  and  jaw  and  chest  injuries. 

3.  The  head  acceleration  curve  is  smoothed  by  the  airbag  with¬ 
out  a  significant  reduction  from  the  peak  values.  However, 
by  limiting  the  excursion  of  the  head,  the  inflatable  collar 
could  be  effective  in  reducing  head  injury  due  to  Impact  on 
the  instrument  panel  or  other  parts  of  a  collapsing  cockpit. 
In  other  words,  the  airbag  does  not  significantly  reduce 
peak  accelerations  in  the  range  up  to  10  g  in  sled 
acceleration.  However,  its  effectiveness  at  higher  accelera¬ 
tions  is  still  an  open  question.  Cadaver  or  animal  tests 
would  be  needed  at  higher  acceleration  levels,  before  risking 
live  subjects  with  inflatable  collars. 

4.  The  "Coulomb  damping"  model  used  for  the  neck  moment  in 
the  analytical  part  of  the  work  shows  good  agreement  with 
experiment  for  acceleration  levels  which  do  not  result  in 
neck  injury  in  human  subjects.  The  solid  damping  in  the 
voluntary  range  of  neck  rotation  plus  nonlinear  hard  springs 
to  stop  the  motion  is  therefore  presented  in  this  report  as 
an  alternative  to  viscoelastic  models. 

5.  The  check  on  the  theory  using  the  experimental  data  reported 
by  Ewing, Thomas,  et.al.,  Ref.  (1),  was  hampered  by  the 
absence  of  plots  of  the  acceleration  of  the  base  of  the  neck 
(T-l  acceleration).  The  reported  sled  acceleration  was 
used  in  checking  out  the  computer  program  for  the  current 
analysis  as  an  approximation  for  the  T-l  acceleration,  but 
this  is  not  satisfactory  for  establishing  realistic  values 
for  the  parameters  involved. 


The  difference  between  the  sled  acceleration  and  the  neck 
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acceleration  is  due  to  the  mechanics  of  the  restraint 
system  of  lap  belt,  shoulder  harness,  and  chest  strap.  This 
difference  was  recognized  by  Ewing,  Thomas,  et.  al.  as 
they  note  the  following:  "The  effect  of  the  restraint 
variable  is  assumed  to  be  independent  of  the  fundamental 
relationship  between  the  T^  acceleration  and  the  head 
acceleration.  However,  as  a  practical  matter,  each 
individual  restraint  system  was  adjusted  as  tightly  as  pos¬ 
sible  for  each  run,  both  as  a  safety  measure  and  to  keep 
the  photographic  targets  within  the  field  of  view  of  the 
sled-mounted  cameras.  Parenthetically,  it  should  be  noted 
that  the  restraint  adjustment  procedure  was  determined 
experimentally  to  result  in  apparent  excellent  control  of  the 
restraint  variable,  as  will  be  shown  later.  An  analysis 
of  the  precise  effects  of  variations  in  the  restraint  of  an 
individual  subject  on  the  dynamic  response  of  the  head  and 
neck  must  be  performed,  but  has  not  yet  been  accomplished." 

Therefore,  the  data  necessary  to  plot  the  T-l  accelerations 
were  recorded  with  the  intent  of  analyzing  it  as  required. 

One  of  the  concluding  remarks  in  Ref.  (1)  is  "However,  the 
variables  presented  in  the  report  are  only  a  partial  sampling 
of  what  is  presently  available.  Exhaustive  and  effective 
use  of  the  data  requires  large  scale  digital  computer  inter¬ 
action  with  the  data  base  on  a  continuing  basis.  The  specific 
interaction  should  be  determined  by  the  specific  problem 
posed." 

It  seems  worth  proposing  that  the  T-l  acceleration  data  be 
plotted  and  compared  with  the  sled  acceleration  data  and 
the  head  acceleration  data. 

6.  The  experimental  run  profile  for  each  subject  in  Ref.  (1) 
contained  runs  with  sled  accelerations  from  3G  to  10G.  The 
plotted  results  are  for  runs  at  levels  greater  than  5.5G. 
Presumably,  this  is  because  the  head  rotation  response  was 
not  as  dramatic  at  the  lower  levels.  However,  this  data 
would  be  useful  in  comparing  the  theory  in  the  present  report 
and  the  usual  viscoelastic  models.  The  solid  damping  model 
would  predict  little  response  below  a  threshhold  level  of 
acceleration  while  a  viscoelastic  model  would  react  to  any 
acceleration. 

7.  The  estimation  of  parameters  in  the  constitutive  relations 
by  using  Newton's  method  provides  a  rational  means  of 
calculation  that  can  be  incorporated  into  digital  computer 
programs  along  with  the  solution  of  the  equations  of  motion. 
The  lack  of  uniqueness  in  the  normal  equations  for  the  least- 
squares  curve  fit  was  a  nuisance,  but  varying  one  constant 

at  a  time  gave  satisfactory  answers.  The  difficulty  in 
inverting  the  normal  equations  for  polynomial  curve  fits 
is  well-known  and  is  one  of  the  reasons  for  using  orthogonal 
polynomials  as  coordinate  functions. 
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On  Che  practical  side,  the  lack  of  uniqueness  is  an  ad¬ 
vantage.  Different  laws  can  be  used  to  fit  the  same  data 
with  good  accuracy.  Therefore,  different  mathematical 
models  can  be  proposed,  e.g.,  the  solid  damping  plus 
nonlinear  springs  compared  to  viscoelastic  models. 

8.  After  the  least-squares  Iteration,  the  computer  program 
calculates  head  and  neck  moments  supplied  by  the  muscles. 

In  the  normal  rotation  range,  these  results  are  comparable 
in  magnitude  to  measured  static  results  and  dynamic  results 
reported  in  the  literature.  This  correlation  supports 
the  assumptions  of  the  theory  derived  here. 

Further  work  is  desirable  in  several  areas: 

1.  The  upper  limits  for  the  present  theory  should  be  determined 
and  an  injury  criterion  developed  for  neck  injuries.  Moment 
and/or  angular  impulse  levels  that  produce  ruptured  discs, 
fractured  vertebrae,  and  torn  ligaments  need  to  be  estab¬ 
lished.  Experimental  data  in  this  area  cannot  be  obtained 
from  human  volunteers.  Therefore,  indirect  methods  using 
cadavers,  animals,  and  mathematical  models  are  necessary. 

2.  Aircraft  structure  accelerations  are  transmitted  to  the  head 
through  the  seat  and  harness  restraint.  The  optimum  design 
of  this  system  can  reduce  peak  head  and  neck  accelerations. 

The  quasilinearization  technique  used  in  the  current  theory  can 
be  extended  to  the  optimization  problem  of  the  system  as 

a  whole. 

3.  lue  trend  in  the  literature  is  to  write  occupant  simulation 
programs  that  predict  gross  motion  of  vehicles  and  passengers 
using  as  few  degrees-of-f reedom  as  possible.  Then,  using  the 
gross  motion  as  input  data,  more  detailed  programs  are  written 
to  describe  selected  parts  of  the  overall  system.  This 
results  in  a  great  saving  in  computer  cost  at  no  loss  in 
physical  understanding  of  the  mechanisms  affecting  the  motion. 
The  computer  program  described  in  this  report  can  be  extended 
to  give  a  more  detailed  simulation  of  the  head  and  neck  and 

it  can  also  be  coupled  to  a  program  for  simulation  of  a 
complete  crash  situation. 

A.  The  experimental  sled  should  be  modified  so  that  the  effect  of 
the  harness  on  the  T-l  impulse  can  be  included.  This  will 
require  that  the  base  of  the  neck  be  allowed  to  move  longi¬ 
tudinally  with  respect  to  the  sled  against  the  retarding 
force  of  a  non-linear  spring.  Slack  in  the  harness  will  also 
be  simulated. 

5.  The  modified  sled  can  be  used  to  study  the  effect  of  harness 
parameters  on  the  motion  of  the  head  under  impulsive  loading. 
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III.  EXPERIMENTAL 


The  objective  of  the  experimental  program,  as  mentioned  earlier, 
was  to  study  the  potential  of  an  inflatable  collar  as  a  means  of  re¬ 
ducing  head  and  neck  rotations  of  naval  flight  personnel  in  the  instance 
of  crash  induced  frontal  accelerations.  Our  previous  program,  Ref.  (19) 
involved  the  development  of  a  small  test  sled  for  imposing  accelerations 
on  a  simulated  full  size  head  and  neck.  This  consisted  of  a  35  lb. 
sled  mounted  on  cylindrical  guides  and  accelerated  by  a  pneumatic 
launching  system  as  shown  in  Figure  1.  A  pneumatic  snubbing  system. 
Figure  2,  brings  the  sled  to  rest.  The  launching  system  provides  a 
square  wave  pulse  up  to  50  g's  in  magnitude  over  a  stroke  of  2  ft. 

The  pneumatic  snubbing  system  produces  a  triangular  stopping  impulse 
over  a  stroke  of  2  ft.  with  adjustable  onset  rate,  peak  acceleration  and 
fall-off  rate.  The  head  may  be  oriented  so  that  the  test  pulse  is 
Imposed  either  in  starting  or  stopping  of  the  sled. 

On  the  previous  program,  a  set  of  experiments  was  conducted  with 
an  Alderson  aluminum  and  rubber  anthropometric  dummy  head  attached  to 
the  sled  with  a  simple  pinned  linkage  the  same  length  as  the  human  neck. 
An  inflated  section  of  fire  hose  was  used  under  the  chin  to  reduce  head 
rotations.  The  simple  hinge  was  an  over  simplification  of  the  neck; 
however,  these  experiments  indicated  that  significant  reductions  in 
head  accelerations  could  be  achieved  if  the  chin  bag  or  an  inflatable 
collar  could  be  developed.  The  test  results  also  provided  a  valuable 
object  of  comparison  for  the  analytical  head  rotation  model  developed 
on  the  program  for  use  in  a  parameter  study. 

The  results  of  the  experimentation  on  the  first  program  indicated 
that  the  current  program  should  address  itself  to  the  development  of  a 
more  realistic  model  of  the  human  neck  so  that  better  experimental 
results  could  be  obtained. 

An  iterative  approach  was  taken  in  the  design  of  the  neck.  After 
consulting  an  anatomy  book,  the  critical  elements  were  simulated  with 
mechanical  parts,  (Figure  3).  The  seven  vertebrae  simulated  with  plastic 
discs,  were  separated  with  discs  of  foam  rubber  and  the  assembly  held 
together  with  str.'inds  of  steel  cable.  The  primary  muscles  were  simula¬ 
ted  with  elastic  cord  covered  with  a  fabric  braid  (called  shock  cord) . 
Four  muscles  were  used;  one  at  each  side  and  two  at  the  back.  Sled 
tests  were  conducted  with  the  head  and  neck  to  evaluate  neck  perfor¬ 
mance.  Accelerometer  and  high  speed  photographic  data  were  compared 
with  data  from  one  of  the  live  subject  tests  from  Ref.  1  (Subject 
004,  6.9  g's)  for  which  a  high-speed  film  was  available.  Adjustments 
were  made  in  the  neck  until  the  motion  of  the  head  and  neck  closely 
simulated  that  of  the  live  subject.  The  final  neck  design  is  shown 
in  Figure  4. 

To  facilitate  the  comparison  of  the  model  test  data  with  live 
subject  data,  a  small  analogue  computer  was  built  to  provide  directly 
the  vector  sum  of  the  radial  and  tangential  head  accelerometers  for 
display  and  recording  on  the  oscilloscope.  This  resulted  in  a  consider¬ 
able  reduction  in  the  effort  required  to  make  the  desired  comparisons. 

The  instrumentation  used  for  the  sled  tests  is  shown  in  Figure  5. 
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A  comparison  of  the  resultant  accelerations  of  the  model  head 
with  that  of  two  live  subjects  is  given  in  Figure  6.  The  neck  base 
impulse  for  the  model  test  was  the  same  as  that  measured  by  an  accelero¬ 
meter  at  the  base  of  the  neck  of  Subject  003,  Ref.  (1)  (Figure  7); 
subject  LMP  of  Ref.  16  was  subjected  to  a  9.6  g  half  sine  wave  impulse.* 
These  curves  have  the  same  characteristic  shape  which,  with  the  exception 
of  the  initial  peak,  closely  resembles  the  shape  of  the  sled  impulse. 

The  time  and  amplitude  of  the  initial  peak  were  found  in  the  model  tests 
to  be  a  function  of  the  head  translation  which  occurs  before  head 
rotation  begins.  A  greater  resistance  to  head  translation  results  in 
less  translation  and  less  jerking  just  before  head  rotation  begins.  It 
was  noticed  that  the  orientation  of  the  neck  base  plate  on  the  model 
influenced  the  size  of  the  initial  peak.  It  is  suggested,  therefore, 
that  the  size  of  the  initial  peak  and  the  time  required  to  achieve  it 
can  be  affected  by  the  posture  of  the  subject  in  the  seat  as  well  as 
the  subject's  general  ability  to  resist  head  translation.  Also,  it 
is  believed  that  the  harness  and  seat  characteristics  have  an  effect 
since  slack  in  this  system  can  serve  to  sharpen  the  pulse  transmitted 
to  the  subject. 

As  previously  indicated,  the  motion  of  the  model  was  compared 
to  the  head  motion  of  subject  004  through  the  use  of  high  speed  film 
sequences  of  the  tests;  such  a  comparison  for  the  final  neck  design  is 
presented  in  Figures  8  and  9.  As  can  be  seen  here,  the  motions  are 
very  similar. 

An  important  feature  of  the  model  neck  is  the  cable  restraint  which 
holds  the  chin  down  when  the  head  rotates  forward.  This  simulates  the 
live  subject's  ability  to  prevent  his  head  from  bobbing.  The  cable 
feeds  freely  through  a  guide  as  the  head  rotates  forward  (Figure  4)  but 
is  secured  by  a  knurled  eccentric  wheel  when  the  direction  of  the  head 
rotation  is  reversed. 

The  effectiveness  of  the  inflatable  collar  concept  was  evaluated 
using  an  inflated  bag  under  the  chin.  This  bag  was  made  from  a  section 
of  fire  hose  as  illustrated  in  Figure  10.  A  series  of  tests  was  made 
with  the  head  facing  the  launcher  so  that  a  testing  impulse  approx¬ 
imating  a  square  wave  was  applied  to  the  base  of  the  neck  during  the 
test  period.  The  resultant  head  acceleration  for  impulses  of  10  g's 
is  plotted  in  Figure  11  for  the  head  without  restraint  and  with  the  chin 
bag  inflated  to  5,  10,  and  25  psi.  A  similar  plot  for  5  g's  is  presented 
in  Figure  12  for  the  unrestrained  head  and  the  head  with  the  25  psi 
chin  bag.  In  general,  without  the  chin  restraint,  the  resultant  head 
acceleration  curve  is  similar  in  shape  to  the  applied  impulse  except 
for  the  initial  peak  and  is  greater  in  magnitude  and  longer  in  duration. 
The  curves  in  Figure  12  indicate  that  the  chin  bag  reduces  the  size  of 


*  Neck  base  accelerometer  data  was  not  available  for  subject  004  so 
the  data  for  subject  003  was  used.  It  is  assumed  that  these  were  similar. 


the  resultant  head  Impulse  but  that  Increases  In  bag  pressure  over  5  psl 
have  little  effect. 

Static  tests  were  also  conducted  to  measure  the  effect  of  the  chin 
restraint.  The  test  set  up  is  illustrated  in  Figure  13.  Weights 
were  used  to  impose  a  horizontal  force  on  the  top  of  the  head  and  head 
rotations  were  measured.  Plots  of  head  rotation  as  a  function  of  applied 
force  are  presented  in  Figure  14.  It  is  interesting  to  note  here  that 
the  greatest  single  effect  at  large  rotations  is  achieved  with  the  5  psi 
inflation  pressure.  While  pressures  up  to  25  psi  were  used  in  the  bag, 
it  is  believed  that  these  higher  bag  pressures  might  pose  a  threat 
to  a  person's  Jaw  structure. 

In  summary,  the  experimental  program  set  out  to  build  an  improved 
model  of  the  human  neck  in  order  to  obtain  data  on  the  potential  of  a 
chin  bag  or  inflatable  collar  as  a  means  of  reducing  head  rotations 
associated  with  crash  induced  frontal  accelerations  on  Navy  flight 
personnel.  A  neck  was  developed  which  closely  simulates  the  behavior 
of  the  human  neck  in  the  desired  mode  of  operation.  This  neck  was 
used  for  experiments  with  an  inflated  chin  bag.  The  chin  bag  test 
results  indicate  that  some  reduction  in  head  accelerations  can  be 
achieved.  Aside  from  chin  bag  tests,  the  mechanical  neck  model  has 
also  been  useful  in  gaining  an  understanding  of  neck  and  head  motions 
and  accelerations  under  impulsive  frontal  loading  similar  to  that 
experienced  by  flight  personnel  in  crash  landings. 
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IV.  THEORY 


The  analytical  work  described  Jn  this  section  of  the  report 
paralleled  the  experimental  work  in  the  previous  section,  but  the 
objective  was  not  to  write  equations  of  motion  for  the  mechanical 
neck.  The  intent  was  to  develop  a  theory  to  describe  the  head  and  neck 
motion  of  humans.  Then,  the  sensitivity  of  parameters  appearing  In  the 
theory  were  evaluated  by  Newton's  method  and  the  parameters  evaluated  by 
a  least-squares  curve  fit  between  the  theory  and  the  experiments  on 
volunteers  reported  by  Ewing  and  Thomas  in  Ref.  1. 

The  emphasis  was  on  the  estimation  theory,  so  the  equations  of 
motions  were  restricted  to  a  two-degree-of-freedom  model  containing 
six  free  parameters.  The  two  degrees  of  freedom  are  head  rotation  and 
neck  rotation. 

The  rigid-body  equations  of  motion  for  the  head  and  neck  are 


(see  Figure  15). 

-F  i  «  mh  a h  (1) 

1 2  +  1.1  ■  %  in  (2) 

-Mlk  +  x  F  l  -  lh  4>  k  (3) 

-M2!l  +  Mlk  +  (£l-P2>  x  Fj-£2  x  F2  -  Ina  k  (4) 

The  kinematic  relations  are 

In  ■  Rn  +  £2  (5> 

Ih  =  Fn  +  £l  +  Ph  (6) 

In  =  in  =  In  +  «  k  x  £2  -  «2£2  <7> 

Ih  =  ih  =  Rn  +  (4x£l  +  ik  x  3  -4>2Ph  (8) 


The  forces  Fj  and  £2  can  be  eliminated  from  Eqs.  (3)  and  (4)  by 
solving  in  terms  of  a^  and  a^  from  Eqs.  (1)  and  (2).  Then  replacing  a^ 
and  ^  results  in  the  differential  equations 

Uh  +  ®hPf!)  $  +  “hPhPi®  cos  (4>-a) 

=  -Mi  +  m^Ph  [xn  cos  $  +  y’n  8in  <J>  “Pl“2  sin  (4>-oi)]  (9) 

"hPhPl^  cos  (<H*>  +  (In  +  mnP2  +  m^)  “ 

*  Mi  -M2  +  mhPhPl  sin  (<M0  +  (n>hPl  +  cos  a  +  yn  sin  cx) 

(10) 

The  acceleration  components  at  the  base  of  the  neck,  xn  and  yn, 
are  given  functions  of  time.  The  mass  properties  and  length  are  anthro¬ 
pometric  data  that  can  be  measured  or  estimated  from  volume  measurements 
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NOTATION  (see  Figure  15) 


Ih 


Moment  of  inertia  of  head  about  its  center  of  mass, 
pound-in-seconds2. 


I 


n 


Moment  of  inertia  of  neck  about  its  center  of  mass, 
pound-in-seconds2 . 


F1-N1  et  -  Qi 
F2“Ni  et  +  Q2 

mi.m2 

m  h*  m  n 

A.  1.  A 

et  «•  -  sin  a  i^ 
e^  -  -  cos  a  i 

a 

Ph 

pl 


ea  Force  between  head  and  neck,  pounds. 

ea  Force  between  neck  and  shoulders,  pounds. 

Moments  at  occipital  condyles  and  moment  at  T^  vertebrae, 
respectively  inch-pounds. 

I 

Mass  of  head  and  neck  respectively,  pounds  -  second^/inch. 

Right-handed  set  of  unit  vectors.  is  vertical  and  k 
is  out-of-plane  of  paper. 

+  cos  a  Unit  vector  fixed  to  neck. 

-  sin  a  Unit  vector  perpendicular  to  neck. 

Angle  of  head  rotation  from  the  vertical  (measured  from  J 
direction),  radians. 

Angle  of  neck  rotation. 

Distance  from  occipital  condyles  to  center  of  mass  of 
the  head,  inches. 

Length  of  neck  from  T^  vertebrae  to  occipital  condyles, 
inches. 


P 


2 


Distance  from  T^  vertebrae  to  center  of  mass  of  the 
neck,  inches. 


Bn 


*n  A  +  i 


Acceleration  of  the  neck  at  Tj  vertebrae,  inches/ 
second 2. 


g 


Acceleration  of  gravity,  inches/second^. 
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and  assumptions  on  mass  density. 

The  key  problem  is  the  const!  tut  ive  relations  that  determine  the 
head  and  neck  moments  and  M2 

Mj  -  Mj  (t,<j),<i>,<j>, «,«,«) 

M2  *  M2  (t,<J>,$,<J>,u,«,u) 

Until  these  relations  are  defined,  the  equations  of  motion  cannot  bo 
integrated  explicitly. 

Several  different  constitutive  relations  have  been  proposed  for 
the  head  and  neck. 

Bowman  and  Robbins  (2)* adopted  the  Maxwell  element  of  a  linear 
spring  and  a  linear  viscous  damper  in  series.  The  spring  and  damper 
coefficients  are  assumed  to  be  linear  functions  of  "the  voluntary 
static  moment".  This  latter  moment  is  not  clearly  defined.  The 
Maxwell  element  was  proposed  by  Moffat,  Harris,  and  Haslam  (3)  based 
on  a  study  of  knee  moments.  They  obtained  coefficients  for  the  knee  by 
impressing  a  sinusoidal  motion  on  the  leg  and  measuring  directly  the 
total  moment  necessary  to  produce  the  motion.  The  resultant  moment 
was  split  into  components  and  the  moment  due  to  the  muscles  was  fit 
to  the  Maxwell  model  using  least-squares.  In  the  experiment,  a 
voluntary  knee  moment  was  measured  for  the  knee  reacting  against  a 
foot  plate  before  the  test  apparatus  was  set  in  motion. 

In  the  notation  of  the  present  report,  the  constitutive  relations 
would  be 


Ml  +  Mi  *  +  ki  (<j>-oi)  (11) 

^2  +  m2  =  +  k2«  (12) 

where  the  kj[  and  cj  coefficients  are  linear  functions  of  initial  muscle 
tension.  This  proposed  relation  has  the  advantage  that  the  angular 
velocity  and  acceleration  resulting  from  an  applied  impulse  do  not 
depend  on  the  initial  values  of  the  angles.  This  independence  from 
initial  values  seems  in  agreement  with  experimental  data.  However, 
the  Maxwell  element  responds  to  any  applied  moment,  no  matter  how  small 
or  slowly  applied.  This  feature  does  not  seem  in  agreement  with  ex¬ 
perimental  observations. 

A  general  linear  viscoelastic  element  similar  to  the  Maxwell 
element  was  used  in  the  analysis  of  whiplash  reported  by  McKenzie 
and  Williams  (4).  The  same  type  of  viscoelastic  elements  are  included 
in  the  analytical  neck  simulation  of  Melvin,  McElhaney,  and  Roberts  (5). 
The  linear  viscoelastic  element  was  first  assumed  by  Orne  and  Liu  (6) 
as  constitutive  relations  for  intervertebral  discs  in  their  mathematical 
model  for  spinal  response  to  axial  impact  simulating  the  pilot  ejection 
problem. 


*  Number  in  parentheses  refer  to  list  of  references. 
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The  linear  viscoelastic  model  would  assume  the  relations 

M1  +  pl^l  “  %1  [ 4>-ct-4> (o)  +  ot(o)]  +  qn  [$-ot]  (13) 

M2  +  P2M2  =  qo2  [a-a(o)]  +  q12a  (14) 

during  the  normal  range  of  motion.  Stops  would  be  included  In  the 
mathematical  model  to  limit  the  angular  motion.  The  linear  spring  con¬ 
stants  qQ^  multiplying  angular  displacement  is  certainly  inappropriate 
for  simulating  the  effect  of  muscles  for  quasi-static  loadings. 

The  constitutive  relations  actually  used  in  the  present  study  were 

M1  =  [^1  +  D1  ($-oi)>0, 

M2  -M1  =  [c2  +  D2an2]  mhg  a>0 

*  -Cimhg  (4>-a)<0 

m2  -Mi  *  [-C2  +  D2an2]  mhg  a<0 

The  equations  are  scaled  by  head  weight,  m^g.  The  equations  are 
designed  to  represent  flexion  such  as  produced  by  the  sled  accelera¬ 
tions  in  the  experiments  described  in  Ref.  (1).  The  constant  terms 
and  C2  oppose  the  angular  motion  no  matter  in  which  direction  the  head 
is  moving,  in  the  same  manner  as  Coulomb  (solid)  friction.  The  non¬ 
linear  springs  are  assumed  to  have  large  positive  exponents,  n^  and 
n2,  so  that  these  hard  spring  terms  become  significant  near  the  limits 
or  voluntary  motion  and  act  effectively  as  stops.  For  hyperextension 
or  whiplash  problems,  Eqs.  (15c)  and  (15d)  would  need  to  be  modified  to 
include  similar  stops  for  negative  values  of  the  head  and  neck  rotation. 
Also,  since  the  neck  is  not  as  strong  in  extension  as  it  is  in  flexion, 
the  constants  Cl  and  C2  should  be  replaced  with  different  values  when 
the  neck  is  extending.  This  point  is  discussed  further  in  the  results 
section  of  this  report. 

As  pointed  out  recently  by  Soechting  and  Paslay  (7),  there  will 
be  a  response  time  before  the  muscles  can  react  to  resist  the  motion. 

The  algebraic  sign  associated  with  the  constant  terms  Cj  and  C2  should 
depend  on  angular  velocities  at  (t-l)  where  T  is  a  neural  delay  time 
rather  than  at  the  current  value  of  time  to  t.  This  delay  time  is  on 
the  order  of  .1  second  for  the  stretch  reflex  of  relaxed  arm  muscles. 

The  reaction  time  of  the  neck  muscles  of  human  subjects  expecting  an 
acceleration  pulse  is  not  known,  however,  if  it  is  on  the  order  of  100 
msec,  this  delay  time  could  be  a  significant  fraction  of  the  total  time 
for  the  event.  One  way  of  measuring  this  effect  on  volunteers  would  be 
to  reproduce  constant  acceleration  levels  with  different  onset  rates. 

A  logical  modification  of  the  constitutive  relations  in  Kqs.  (15) 
is  to  include  a  delay  time.  However,  the  results  presented  later  do  not 
contain  this  effect. 


(4>— nt)  >0  (15a) 

(15b) 

(15c) 

(15d) 


Aside  from  differences  in  the  assumed  form  of  the  constitutive 


relations  in  biomechanics,  there  is  general  agreement  that  experimental 
data  are  lacking  to  fix  the  values  of  the  constant  parameters  that  appear 
in  the  equations.  "Appropriate  parameter  values  for  the  neck  are  not 
available  In  the  literature,  but  Moffat,  Harris  and  Has  lam  (9)  give 
experimentally  determined  values  for  the  knee  Joint.  These  values  were 
used  following  application  of  appropriate  scaling  procedures  using  other 
available  experimental  data  to  obtain  rough,  order-or-magnitude  values 
for  corresponding  neck  quantities."* 

"In  particular,  it  should  be  noted  that  the  values  of  the  material 
properties  governing  the  viscoelastic  behavior  of  the  intervertebral 
discs  were  those  used  by  Orne  and  Liu  (1971).  These  authors  found  that 
a  wide  range  and  combination  of  material  parameters  had  little  effect  on 
the  magnitude  of  the  peak  force  response."** 

"...  muscle  and  neural  feedback  parameters  are  difficult  to 
obtain.  The  time  constant  D  appears  to  be  reasonably  well  established, 
however,  values  for  the  muscle  spring  constant  K  and  the  feedback  para¬ 
meter  C  have  been  determined  only  for  Isolated  sets  of  muscles."*** 

In  addition  to  the  lack  of  data,  the  usefulness  of  certain  basic- 
studies,  e.g..  Ref.  (3)  and  Ref.  (10),  has  been  impaired  by  the  authors 
fitting  the  results  to  a  mathematical  model  rather  than  presenting  the 
raw  data.  There  is  no  objection  to  trying  to  analyze  data,  but  it  should 
not  be  at  the  expense  of  obscuring  the  basic  measurements. 

Since  there  is  a  lack  of  data  on  proper  parameters  for  constitu¬ 
tive  relations  for  muscles  and  joints,  it  seems  necessary  to  adopt  an 
indirect  approach.  This  approach  consists  of  formulating  the  theoretical 
model  and  then  adjusting  the  parameters  to  fit  the  data  of  tests  such 
as  Ref.  (1).  If  the  theory  is  consistent,  the  known  parameters  can  be 
used  to  predict  responses  to  sets  of  applied  forces  different  from  the 
test  data.  This  approach  is  commonly  used  in  mechanics;  in  a  particular 
example  in  biomechanics,  Melvin,  McF.lhaney,  and  Roberts,  Ref.  (5),  fit 
their  mathematical  model  to  the  tests  of  Mertz  and  Patrick  (11). 

The  main  aim  of  the  present  study  has  been  to  develop  a  systematic 
method  for  arriving  at  the  parameters  in  the  constitutive  equations  so 
that  the  theory  fits  the  experimental  data  of  Ref.  (1).  The  result  has 
been  a  computer  program  that  automates  the  search  for  parameters  that  fit 
the  data.  The  program  is  an  improvement  over  varying  different  para¬ 
meters  in  a  random  fashion  and  comparing  the  resulting  motion  to  the 
experimental  results.  As  will  be  shown,  the  current  version  of  the  program 
is  probably  not  an  optimum  algorithm  in  terms  of  computer  time,  and  in 
some  cases,  the  search  procedure  does  not  converge.  However,  the  progress 
thus  far  has  shed  important  light  on  the  interpretation  of  the  experi¬ 
mental  results  of  Ref.  (1)  and  could  serve  as  a  guide  in  conducting 


*  Ref.  (2)  Bowman  and  Robbins 

**  Ref.  (4)  McKenzie  and  Williams 

***  Ref.  (7)  Soechting  and  Paslay 


future  experiments.  The  theoretical  results  to  date  have  also  suggested 
improvements  In  the  analytical  model. 

The  estimation  of  parameters  in  the  constitutive  relations  is 
treated  by  finding  the  effect  of  varying  each  parameter  by  Newton's 
method.  Ref.  (12)  and  Ref.  (13).  The  variation  in  eacli  parameter  is 
computed  by  a  least-squares  fit  of  the  experimental  data.  This  method 
of  estimation  theory  has  been  used  for  process  dynamics  in  chemical 
engineering  by  Lee,  Ref.  (14),  and  others.  Bellman,  Ref.  (15),  describes 
other  applications,  but  the  method  does  not  seem  to  be  in  common  use  for 
the  biomechanics  of  human  motion. 


The  procedure  is  an  iterative  one.  First,  a  trial  set  of  para¬ 
meters  is  assumed  along  with  the  initial  conditions  and  the  rest  of  the 
data  to  solve  the  nonlinear  equations  of  motion,  Eqs.  (9)  and  (10),  for 
<p  -  4>0(t)  and  a*ao(t).  This  initial  solution  is  calculated  in  the  computer 
program  by  a  Runge-Kutta  subroutine.  The  next  step  is  to  compute  the 
effect  of  a  change  in  each  parameter  by  means  of  a  Taylor  series  expan¬ 
sion  about  the  initial  solution.  <J>G(t)  and  ct0(t). 

<J>  *  d>o  (16a) 

(**  %  +  6a  (16b) 


Before  the  Integration  of  the  initial  solution  was  programmed,  the 
equations  of  motion  were  normalized  by  dividing  through  by  the  mass  of 
the  head,  m^.  The  constitutive  relations  are  scaled  by  dividing  the 
moments  by  the  weight  of  the  head,  m^g,  as  indicated  in  Eqs.  (15).  The 
expansion  in  the  parameters  is  then  of  the  form 

Mj  Mi(o)6Mi  (17) 

mh  =  mh  +  mh 


where 


etc. 


(4.0-ao)  >0, 


(<J>-a)>o 


6M1  / 

+  6D1  (<t>o-(*o)ni  +  n1D1(4>o-Oo)nl“1  (6«|>-6a) 

+  6n!  [ln(<J>0-a0)]  } 


Note  that  the  expansion  is  truncated  at  first  order  terms  in  the  cor¬ 
rections.  For  example,  the  term  Di(<J)-a)nl  is  expanded  (d|°^  +  6D^) 

(4>0+6(J)-a0-6a)nl+6nl  =  (4>0-a0)nl  +  6Di(4>0-a0)ni  +  Di(<p0-ao)nl-1 

(6<j>-6a)  +  6ni  ln(4>0-ao)  (<)>0-a0)nl  +  .  .  .  and  quadratic  and  higher  terms 
in  the  corrections  and  their  products  are  dropped.  Similarly,  in 
substituting  Eqs.  (16)  into  the  equations  of  motion,  only  linear  terms 


» 
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are  retained,  such  as 

4>2  sin  (<t>-u)  =  02  sin  (<f>0-«0)  +  24>06<t>  sin  (4>o-oto)  +  cos  ^o-nto^  (<S<l»-«Stx) 

The  resulting  "variational"  equations  of  motion  can  be  written  as 

Lu(6<j>)  +  L12  (6a)  =  I  g(<J>  ,a  )6Ci  (18a) 

AA  1=1  1J  o  o 

L2i(64>)  +  L22  (6a)  =  E  R2i(^0*ao^Ci  (18b) 

where  the  six  subscripted  variables  6Cj  are  6C1 ,  6C2,  6D^,  6D2,6ni,  and 
6n2  in  that  order.  The  Ly  are  second  order  linear  differential 
operators. 

Since  Eqs.  (18)  are  linear  differential  equations,  a  particular 
solution  (6i}>i,6ai)  can  be  obtained  for  each  term  in  the  series  on  the 
right  and  the  results  superposed  to  obtain 

6 

4>  =  4>o  +6<t>  =  4-0  +  Zm  6Cj  64>i  (19a) 

ft 

a  =  a  +  6a  =  aQ  +  E  ^6Cj  6at  (19b) 

The  correction  terms  6 4>  and  6a,  are  linear  functions  of  the  6C1 
so  that  the  original  Intent  was  to  use  Eqs.  (19)  in  a  linear  regression 
function  and  to  compute  the  6C1  by  the  method  of  least-squares  using  the 
experimental  data  of  Ref.  (1). 

Angular  accelerations  would  seem  intuitively  to  be  the  choice 
for  the  linear  regression  function  since  the  moment  terms  affect  the 
angular  accelerations  directly.  However,  there  were  problems  inter¬ 
preting  the  acceleration  data.  These  problems  will  be  discussed  in 
more  detail  later. 

The  regression  function  was  programmed  to  fit  the  measured 
angular  deflections.  The  idea  was  to  solve  the  "normal  equations"  for 
the  six  6Ci,  solve  Eqs.  (19)  explicitly  for  4>  and  a,  repeat  the 
iteration  with  4>  and  a  replac ing  4>o  and  aQ  respectively  in  Eqs.  (16). 

The  iteration  continues  until  the  corrections  6Cj  become  small.  However, 
this  procedure  failed  to  converge  because  the  matrix  of  the  normal 
equations  was  nearly  singular.  This  difficulty  took  considerable  time 
to  overcome  because  it  is  not  clear  what  caused  the  trouble.  Are  the 
coordinate  functions,  64>j  and  6aj  in  Eqs.  (19),  linearly  dependent? 

Does  the  numerical  analysis  contain  truncation  and  round-off  errors 
that  affect  the  inversion  of  the  normal  equations?  Is  the  difficulty 
a  combination  of  using  nearly  linearly  dependent  functions  and  numerical 
errors? 

Hamming,  R.  and  Feigenbaum  (17)  explain  the  situation  by  saying 
"It  frequently  happens  that  the  system  of  equations  does  not  determine 
the  parameters  very  accurately.  This  tends  to  worry  the  beginner.  The 
reason  they  are  poorly  determined  is  often  that  the  minimum  is  broad 
and  rather  flat.  It  is  true  that  the  optimum  values  are  then  poorly 
known,  but  it  is  also  true  that  whichever  set  you  choose  (among  those 
that  look  good)  does  not  very  much  affect  the  value  of  the  function  you 
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were  minimizing;  thus,  the  uncertainty  does  you  little  harm." 

As  an  example,  they  cite  the  1  Inear  regression  problem  using 
powers  of  x.  They  are  linearly  independent  functions;  but  as  the  degree 
of  the  regression  polynomial  get  larger,  the  accurate  inversion  of  the 
normal  equations  becomes  more  difficult. 

It  Is  not  difficult  to  construct  counter  examples  where  the  normal 
equations  are  actually  singular.  The  simple  linear  oscillator  is  a 
familiar  case  that  illustrates  this  point.  The  steps  in  the  numerical 
solution  can  be  shown  explicitly  so  it  makes  a  good  illustrative  example. 
The  differential  equation 

mx  +  kx  *  0  (20) 

governs  the  motion.  Assume  that  the  mass  m  is  given  and  an  experimental 
curve 


xe  ■  A  sin  uit 

is  available.  The  problem  is  to  find  the  spring  constant  k  by  Newton's 
method.  In  general,  the  linear  variational  equations  have  variable 
coefficients  and  can't  be  solved  in  closed  form,  but  can  be  solved 
numerically. 

The  procedure  is  to  assume  an  initial  value  of  k  =  kQ  and  seek  a 
correction  6k  so  that  k  =  k0  +  6k  gives  an  improved  approximation  to  the  ex¬ 
perimental  data.  The  solution  of  Eq.  (20)  for  k  *  k<>  is  denoted  by  x*x0. 

Whe;  k  is  varied  by  6k  the  variation  in  x  is  6x.  Substituting  in  Eq.(20)  gives 

m(xQ  +  6x)  +  (kQ  +  6k) (x0  +  6x)  *  0  (21) 

Dropping  the  nonlinear  term  in  the  variations  and  using  the  fact  that 


by  definition 

mx'o  +  k0x0  *  0  (22a) 

2 

x0  *  AQ  sin  0Jpt  0)o  *  kD/m  (22b) 

leads  to  the  variational  equation 

m6x  +  k06x  »  -6k  Ap  sin  uipt  (23) 

The  solution  of  Eq.  (23)  is 

6x  *  — —  6k  t  cos  ojQt  (24) 

2oj0m 

The  current  approximation  for  x  is 


x  “  xQ  +  6x  =  A0  sin  uiQt  +  Ji a.. .  6k  t  cos  a>0 1 

2(Uom 


This  approximation  is  compared  to  the  experimental  curve  xe  by  a  linear 
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regression  analysis  to  determine  the  correction  6k 


Ai- 

2(u0m 


sin  u)t 


-A0  sin  di0t) 


o  o 

t  cos  u)0t  dt  (25) 

The  integrals  in  Eq.  (25  )  depend  on  the  upper  limit  T,  but  the  leading 
term  for  6k  does  not. 

6k  =  2<jJ0m(u)-<d0)  +  .  .  .  (26) 


Setting 


kj  =  kQ  +  6k 


the  iteration  can  be  repeated.  In  general,  equations  such  as  Eqs.  (22a), 
Eq.  (23)  and  Eq.  (25)  must  be  integrated  numerically.  Tn  this  simple 
example,  Eq.  (26)  for  the  correction  on  kG  can  be  checked  directly  since 
the  solution  of  the  original  problem  can  be  written  in  closed  form. 


ur 


k  _  kp  +  6k  _  2  +  -- 

ra  m  o  T  tn 


2 

6k  =  m(urt-w0)  (w-«0)  =  2mt<)0(<i>-ui0)  +  m(aj-U)e)  (27) 

Comparing  Eqs.  (26)  and  Eqs.  (27)  shows  that  the  approximate  expression 
differs  from  the  exact  expression  by  a  quadratic  term  in  the  error  in 
the  current  approximation  to  the  angular  frequency  so  that  convergence 
is  rapid  and  the  correct  value  of  k  is  determined  easily.  However,  if 
the  mass  m  is  unknown  as  well  as  the  spring  constant  k,  the  procedure 
breaks  down.  Varying  both  m  and  k  results  in  an  expression  for  the 
variation  in  x  as 

6x  =  — ^ —  6k  t  cos  wt  -  Aoafe  6m  t  cos  u>t 
2a)omo  2m0 

The  pair  of  linear  normal  equations  equivalent  to  Eq.  (25)  are  obviously 
singular.  In  the  closed  form  solution,  it  can  be  seen  that  the  ratio 
k/m  is  a  single  independent  parameter  and  the  singular  normal  equations 
reflect  this. 

The  mass  of  the  head  m^  is  used  to  normalize  Eqs.  (9),  (10),  and 
(15)  so  that  this  particular  variable  is  not  redundant  in  the  analysis 
of  head  and  neck  motion.  In  the  example  of  the  linear  oscillator,  the 
rank  of  the  matrix  of  the  normal  equations  is  one  or  one  less  than  the 
number  of  variables  when  the  mass  and  spring  constant  are  both  varied. 
During  the  checkout  of  the  computer  program  for  head  and  neck  motion, 
different  sub-matrices  of  the  normal  equations  were  inverted  to  see  if 
sub-sets  of  the  six  parameters  could  be  varied  independently.  In  some 
cases,  they  could  and  in  others,  the  sub-matrices  were  also  nearly 
singular.  The  final  version  of  the  computer  program  varies  only  one 
parameter  at  a  time,  however,  the  possibility  of  varying  more  than  one 
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but  less  than  six  In  a  systematic  fashion  Is  a  possible  variation  of  the 
program. 

Before  leaving  the  example  of  the  linear  oscillator,  another 
feature  of  that  solution  should  be  mentioned.  The  particular  solution, 
Eq.  (24),  of  the  variational  equation  grows  linearly  with  time.  This 
linear  growth  affects  the  convergence  of  the  normal  equations,  Eq.  (25), 
if  the  interval  of  integration  T  is  larger  than  several  cycles.  This 
nonuniform  convergence  can  be  avoided  if  a  weighting  factor  t-^  is 
introduced  in  the  error  analysis.  A  weighting  factor  of  t”l  was  intro¬ 
duced  arbitrarily  in  the  regression  analysis  for  the  head  and  neck  but 
it  did  not  improve  the  convergence  and  was  subsequently  removed.  All 
variational  equations  do  not  have  particular  solutions  with  amplitudes 
that  grow  linearly  with  time  so  that  this  feature  of  autonomous  systems 
with  periodic  solutions  is  not  always  present,  Ref.  (18). 

Another  possible  source  of  error  was  checked.  The  six  solutions 
of  Eqs.  (18)  are  assumed  to  be  particular  solutions.  However,  since  the 
differential  equations  are  solved  numerically,  it  is  possible  to  pick  up 
an  unstable  homogeneous  solution  as  part  of  the  numerical  solution.  If 
the  homogeneous  solution  is  dominant,  this  will  affect  the  linear 
independence  of  tl\e  six  solutions  used  in  the  regression  analysis.  A 
seventh  homogeneous  solution  of  Eqs.  (18)  was  computed  numerically  for 
nonzero  initial  conditions.  The  six  particular  solutions  were  made 
orthogonal  to  the  homogeneous  solution.  The  resulting  six  solutions 
were  essentially  unchanged  so  that  this  check  was  also  removed  from 
the  program. 

The  conclusion  from  the  checks  run  on  the  program  is  that  the 
results  are  insensitive  to  the  six  parameters  in  the  sense  that  more  than 
one  combination  of  parameters  will  give  virtually  the  same  response 
curve.  This  conclusion  was  also  reached  by  Orne  and  Liu,  Ref.  (6),  for 
their  viscoelastic  model  of  the  spine.  The  current  version  of  the 
computer  program  varies  one  parameter  at  a  time.  The  parameter  that  is 
varied  is  computed  in  the  analysis.  A  regression  analysis  is  run  for 
the  effect  of  the  variation  of  each  parameter  while  the  others  are 
held  constant  and  the  fractional  change  of  each  parameter  computed.  The 
parameter  with  the  smallest  change  based  on  the  linear  analysis  is 
selected  to  be  changed.  If  the  fractional  change  is  less  than  a 
small  input  number,  the  parameter  with  the  next  smallest  fractional 
change  is  varied.  The  iteration  is  continued  until  all  the  changes  are 
less  than  the  tolerance  number  or  until  an  iteration  counter  is 
exceeded.  The  results  are  then  printed  and  plotted. 
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V.  RESULTS 


Typical  results  are  plotted  in  Figures  16-45  .  Each  computer  run 
generates  five  graphs.  The  runs  are  identified  by  the  subject  number 
and  peak  acceleration  from  Ref.  (1)  and  by  the  six  parameters,  Cj ,  l)j  , 
nj_,  C 2 ,  I>2  and  n£,  that  appear  In  the  constitutive  relations,  Eqs.  (15). 

The  first  two  plots  show  the  head  rotation  and  neck  rotation  as  a 
function  of  time.  These  are  the  two  variables  in  the  regression  analysis 
that  are  fitted  to  the  experimental  data  of  Ref.  (1).  The  experimental 
points  used  as  input  to  the  program  are  also  plotted.  The  experimental 
data  in  Ref.  (1)  extends  over  a  longer  time  period  than  the  points  shown 
on  the  plots.  It  was  found  that  the  head  rotation  reaches  a  maximum 
and  then  rebounds  to  a  more  or  less  constant  value.  The  numerical 
analysis  consistently  arrives  at  the  equilibrium  condition  at  a  higher 
angle  than  experiment.  There  are  at  least  two  possible  explanations 
for  this.  First,  it  is  known  that  the  neck  is  stronger  in  flexion 
than  in  extension.  Ref.  (11).  Therefore,  the  constants  in  Eqs.  (15c) 
and  Eqs.  (15d)  should  be  modified  for  flexion.  This  would  essentially 
introduce  additional  parameters  which  could  be  used  to  fit  the  data  on 
the  rebound  part  of  the  curve.  A  more  direct  way  of  arriving  at  these 
constants  would  be  with  hyperextension  tests. 

The  second  explanation  of  the  apparent  higher  strength  of  the 
model  In  extension  is  that  a  delay  time  should  be  introduced  in  the 
constitutive  relations,  Eqs.  (15).  The  subject  must  sense  the  rebound 
and  tense  opposing  muscles  to  resist  it.  Meanwhile,  the  head  continues 
to  rebound  until  the  subject  can  bring  it  to  rest. 

Whichever  expianation  is  valid,  the  present  theory  is  restricted 
to  flexion  in  one  direction  and  the  curve  fitting  was  done  on  this 
part  of  the  curve  where  the  maximum  velocities  and  accelerations  occur. 

The  head  moment  about  the  occipital  condyles  and  the  neck  moment 
about  Tj  are  the  third  and  fourth  plots.  These  are  normalized  by 
division  by  the  head  weight  so  that  the  units  on  the  plots  are  inch- 
pounds  per  pound.  An  average  head  weight  is  about  10  pounds;  to  convert 
the  plotted  values  for  this  weight  to  foot-pounds,  multiply  by  10/12. 

These  computed  values  of  the  moments  shown  in  the  plots  are  in  good 
agreement  with  measured  dynamic  and  static  moments  reported  in  the 
literature. 

The  last  graph  is  the  computed  resultant  acceleration  of  the  center 
of  gravity  of  the  head.  This  plot  shows  a  peak  acceleration  on  the  order 
of  twice  the  peak  sled  acceleration. 

It  was  intended  to  do  the  regression  analysis  on  the  measured 
head  accelerations,  but  the  experimental  data  from  Ref.  (1)  show  two 
peaks  in  the  acceleration  curve.  The  two  relative  maxima  are  approximately 
the  same  magnitude.  A  typical  experimental  plot  is  shown  in  Fig.  6  of 
this  report. 

The  first  acceleration  peak  is  not  as  pronounced  in  the  analytical 
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results  and  is  of  much  lower  magnitude  than  the  experimental  result. 

The  agreement  between  the  theory  and  experiment  is  better  for  the  second 
peak.  The  computed  amplitude  is  consistently  lower,  but  the  discrepancy 
between  theory  and  experiment  is  less  than  for  the  first  peak. 

As  mentioned  earlier,  the  analysis  uses  the  sled  acceleration  as 
input  for  the  neck  acceleration  since  the  neck  acceleration  data  were 
not  readily  available.  It  seems  likely  that  the  peak  acceleration  at 
the  aeck  of  the  experimental  subjects  was  considerably  higher  than  the 
peax  sied  acceleration  and  occurs  later  in  time.  The  elasticity  and 
slack  in  the  shoulder  harness  and  torso  which  transmit  the  acceleration 
from  the  sled  to  the  upper  torso  act  to  prevent  the  neck  from  having  an 
identical  acceleration  with  the  sled. 

The  one  figure  of  raw  data  in  Ref.  (1),  Figure  DP-4,,  tends  to 
support  this  assumption.  Shulman,  et.al..  Ref.  (20),  measured  a  time 
lag  between  the  beginning  sled  and  subject  accelerations  for  negative 
vertical  acceleration  of  a  seated  pilot. 

The  analytical  model  does  not  have  an  axial  degree-of-freedom. 

The  slow  motion  film  of  one  of  the  tests  reported  in  Ref.  (1)  shows 
an  apparent  stretching  of  the  neck.  An  acceleration  in  this  mode 
could  be  an  alternate  explanation  of  the  first  acceleration  peak 
exhibited  by  the  experimental  data.  However,  the  computed  HtretchinR  of 
the  neck  plotted  In  Ref.  (1)  does  not  indicate  that  this  is  a  signifi¬ 
cant  variable. 

Whatever  the  explanation  for  the  first  acceleration  peak  in  the 
experimental  data,  it  would  be  desirable  to  have  the  neck  acceleration 
data  to  use  as  input  data  for  the  computer  program  rather  than  the  sled 
accelration. 

The  plots  in  this  report  are  based  on  data  for  subjects  3  and  10 
from  Ref.  (1).  For  subject  3,  the  constants  in  the  constitutive  rela¬ 
tion  were  computed  by  the  computer  program  for  each  sled  acceleration 
level.  For  subject  10,  the  constants  were  computed  by  the  regres¬ 
sion  analysis  for  the  7.1gsled  acceleration.  As  a  check  on  the  theory, 
the  other  runs  on  subject  10  were  made  with  the  same  constants. 

The  poorest  agreement  between  theory  and  experiment  is  for  subject 
3  at  9.5  g  acceleratic  n.  It  is  not  apparent  whether  different  initial 
choices  of  the  six  constants  in  the  constitutive  relations  would  improve 
the  regression  analysis  or  whether  the  limitations  in  the  input  data 
and  the  theory  that  have  been  discussed  earlier  are  more  pronounced  at 
higher  acceleration  levels. 

In  summary,  the  results  thus  far  are  encouraging.  They  support 
the  assumption  of  this  report  and  of  other  investigations  that  the  neck 
muscles  affect  the  head  motion  at  survivable  acceleration  levels.  A 
"solid  friction"  type  of  constitutive  relation  for  the  moments  fits 
the  experimental  data  as  well  as  the  viscoelastic  models  found  in  the 
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literature,  although  a  delay  time  should  be  added  to  the  theory  to  allow 
for  the  reaction  time  of  the  subject. 
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VII. 


APPENDIX 


Certain  anthropometric  data  is  necessary  as  input  data  for  the 
computer  program.  These  numbers  were  not  readily  available  for  the 
subjects  of  Ref.  (1).  One  set  of  assumed  data  was  used  for  all  the 
calculations.  The  numerical  values  are  the  following: 

(Ih/mh)  +  pH  •  16.18  m2,  (In+n^p2)  /  mh  -  2.66  in2 

Ph  *  3.0  in..  Pi  *  4.5  in.,  P2  ■  2.25  in., 

mjj/mj,  =  .318 
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FIGURE  5 

Instrumentation  schematic 
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FIGURE  17 


HEAD  MOMENT  PER  HEAD  WEIGHT, IN-LB  PER  LB 
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FIGURE  18 
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SULTRNT  HEAD  ACCELERTflTION.  G'S 
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FIGURE  33 
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FIGURE  44 
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